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Abstract 

1 tie probability density functions (pdf's) of the fluctuating velocity components, 
as well as their first and second derivatives, are calculated using data from the direct 
numerical simulations (DNS) of fully developed turbulent channel flow. It is observed 
that, beyond the buffer region, the pdf of each of these quantities is independent of 
the distance from the channel wall. It is further observed that, beyond the buffer 
region, the pdf s for all the first derivatives collapse onto a single universal curve and 
those of the second derivatives also collapse onto another universal curve, irrespective 
of the distance from the wall. 1 he kinetic-energy dissipation rate exhibits log normal 
behaviour. 


1 Introduction 


It has long been known that the probability density functions (pdf’s) of the velocity fluctua- 
tions in isotiopic turbulence are Gaussian, whereas those of the velocity derivatives deviate 
significantly from Gaussian. This behaviour in grid generated turbulence was first observed 
experimentally by Batchelor and Townsend (1949), and later on by Kuo and Corrsin (1971). 
This fact was later confirmed by She, Jackson and Orszag (1988) using the data from direct 
numerical simulations of isotropic turbulence. The first derivatives of the velocities are ob- 
served to have exponential tails, Narasimha (1989). However, the dissipation rate of kinetic 
energy, which is the sum of squares of the velocity derivatives in different directions, is known 
to exhibit log normal behaviour. 'This was verified experimentally by Kuo and Corrsin (1971 ) 
and by several other investigators referred in their work. 

'Hie presence of exponential tails in the pdf’s of the derivatives implies greater intermit- 
tcncy in the small scale structures. The fact that, even though the pdf’s of the velocities are 



close to Gaussian, the derivative pdf’s have exponential tails, implies that the behaviour of 
the small scale structures is not influenced significantly by the presence of the larger scale 
structures. The knowledge about the pdf of dissipation rate gives us a better picture about 
the intermittency of small scales and may help in developing or assessing the turbulence 
models. The modeled moments for the dissipation rate can be compared with those of the 
log normal pdf. Ascertaining the shapes of the pdf’s of the velocities and their derivatives 
is expected, at a later stage, to help in the development and assessment of the stochastic 
turbulence models. 

The aim of the present work is to observe the nature of the probability distributions of 
the velocity fluctuations, their derivatives and the energy dissipation rate in wall-bounded 
turbulence. In particular, we would like to examine whether the pdf’s depend on the distance 
from the wall or whether they are universal with respect to the wall normal direction. Such a 
universality of the pdf’s was shown by Balachandai and Sirovich (1991) in Rayleigh- Benard 

convection problem. 

In Section 2, a brief description of the DNS data, the coordinate system and the pdf’s 
is given. In Section 3, the question of independence with respect to wall normal direction 
of the velocity pdf’s is examined. In Section 4, the pdf’s of the first derivatives of the 
velocities are examined. In particular the pdf’s are examined for independence of wall 
normal direction and the behaviour of the tails. In Section 5, the second derivatives are 
examined for independence in the wall normal direction. In Section 6, the log normality of 
the dissipation rate is examined. We conclude with a brief summary in Section 7. 

2 Preliminary Considerations 

In our discussions we follow the conventional coordinate system in which x is the streamwise 
direction, y is the wall normal direction and z is the spanwise direction. We will also use the 
indices 1, 2 and 3 to denote .r, y and 7 directions, respectively, and a comma (,) followed by 
a subscript denotes a derivative with respect to the subscript. 
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1 he data used in this study are from temporal DNS of turbulent incompressible channel 
flow. In temporal simulations the flow is assumed to be periodic in the streamwise and 
spamvise directions and to evolve in time. Hence the streamwise and spanwise directions are 
the homogeneous directions. The DNS data of Zang at a Reynolds number of 318 based on 
the channel half height and the wall friction velocity, used in Dinavahi and Zang (1991), are 
used m the present study. Five y locations arc selected for pdf calculations. They are selected 
so as to get one point in the viscous sublayer (y < ?/.,), one in the buffer region between the 
sublayer and the log layer (y < y s ), two in the log layer, and one on the center line. The y 
locations can be selected easily from a plot of a+ vs y+ and are shown as symbols in Figure 
1 , where u + = vju T , i/ + = yu T /v and n r , and v are the wall friction velocity and kinematic 
viscosity, respectively. 1 he value of y s is such that for any y > y s the point is outside the 
viscous sublayer and the buffer region (see Fig. 1). The y+ values in the increasing order 
are: 3.84, 10.12, 63.32, 100.10, and 318.0. By symmetry each selected y corresponds to two 
horizontal planes: one measured from the bottom plate and one from the top plate, thus 
doubling the number of points. There are 216 points in the streamwise direction, 108 points 
in the spanwise direction and 162 points in the wall normal direction. Spanwise symmetry is 
exploited in these calculations. Twenty realizations from the fully turbulent region are used 
in the present calculations. We thus have a total of 20 x ( 2 x 216 x 108) = 933,120 points 
at a given y location to calculate pdf’s. These are close to a million points and are more 
than sufficient to give well resolved pdf’s. 

In order to verify the conclusions arrived at by using the above data, the DNS data of 
Kim (1989) at a Reynolds number of 395 are also used. There are 256 and 192 points in the 
streamwise and spanwise directions, respectively, and 192 points in the wall normal direction. 
No symmetry in the spanwise direction is used in these calculations. Only one realization is 
used in our study and so there are 2 x 256 x 192 = 98,304 points at each y location. The 

following y+ locations are selected: 6.24, 20.50, 122.21, 204.17, 386.24. These locations are 
shown in Figure 2 as symbols. 
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The pdf’s are normalized by the local mean and standard deviation at a given y location. 

In the current simulations, x and 2 being the homogeneous directions, the pdf of a random 
variable <j> is only a function of the wall normal direction and the Reynolds number of the 
flow; hence, it can be denoted as /,(*; », /*„), where <j> is the random variable in the physical 
space, i is the corresponding sample space variable, and R e is the Reynolds number of the 
flow under consideration. This indicates the dependence of pdf on the physical space variable 

y and the parameter R e - 

3 Velocities 

In this section, the pdf’s of the fluctuating velocities are examined as a function of the 
wall normal distance. In all the pdf plots the Gaussian curve has also been plotted as a 
reference for comparison with the calculated pdf’s and the abscissa V is the random variable 
normalized by the local mean and standard deviation. In Figure 3 the pdf’s of the u velocity, 
f u (u\ y,R e = 313), are plotted for three y-locations outside the viscous sublayer. This plot 
shows that the pdf for the u velocity is independent of the y location. These pdf’s are slightly 
skewed to the left. This implies that at a given Reynolds number a single pdf / U (u; R e ) may 
be constructed for points outside the buffer region. From Figures 4 and 5 we come to the 
conclusion that the pdf’s for u and w, f v (v) and f w {w), are also independent of y location 

(for y > y 3 ). 

The pdf’s of «, v and re for Reynolds number 395 are plotted in Figures 6, 7 and 8, 
respectively, and it can be seen that these pdf’s are also independent of y. 

Hence, it seems that the pdf’s of u, e and re, for the two Reynolds numbers considered, 
are independent of y (y > y,) locations beyond the buffer region, i.e., f Ui {u» V, Re) = /«,(«<) 
for y > y 3 - Also, the pdfs show only minor deviations from the Gaussian curve, the slight 

left skewness of the pdf for u being the largest. 


4 


4 First Derivatives of the Velocities 


In this section, we examine the question of independence of y location for the pdf’s of the 
first derivatives. We jwill also consider the hypothesis by Narasirnha. (1989) that the pdf’s of 
the first derivatives are made up of Gaussian core and exponential tails. 

In Figure 9, pdf’s for t/, y at three y locations outside the sub layer are plotted from the 
data at a Reynolds number of 318. From this figure it is seen that the pdf’s of u, y collapse 
onto one single curve for y > y s . Similarity, the pdf’s of the other first derivatives are also 
independent of y location. Having established the independence of y location for the first 
derivatives, pdf’s are calculated by using the entire data set from the three y locations for 
each of the nine derivatives u,*j ( i, j— 1, 2, 3). These pdf’s can be denoted as / U| (5,-^; /? e ), 
showing the independence of y for y > y s , and arc plotted in Figure 10. It can be seen 
that the pdf’s for all nine derivatives fall on the same curve, albeit with some spread in the 
tails. Upon closer examination, it looks like the nine curves may be divided into two groups 
depending on their proximity to each other in tails. The quantities, u, T , u, y , and tc, ? , that go 
into the continuity equation group together and the other derivatives which go into defining 
the vorticity components make another group. This will be investigated in more detail in 
the future. 

To sec if these observations also hold at another Reynolds number, pdf’s for the quantity 
v, x are plotted in Figure 11, using the data at R e = 395. It can be observed that the pdf’s 
beyond the buffer region collapse onto one curve. Analogous to Figure 10 at R e = 318, pdf’s of 
all the nine first derivatives at /i e — 395 are plotted in Figure 12. Figures 10 and 12 suggests 
that the pdf’s of the first derivatives are not only independent of the y location, but also, in 
fact, all of them even fall onto the same curve at a given Reynolds number. 

Having established the independence of y and the universality of the first derivatives, a 
least squares fit to the function 

f{x) = ae- x2/ ° 2 + be-°\ T \ 
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( 1 ) 



is determined for the pdf’s of all the nine derivatives at each of the two Reynolds numbers. 
The above shape is that of a Gaussian core and exponential tails, as suggested by Narasimha 
(1989). The two calculated curves and the experimental data of Van Atta and Chen (1970) 
are plotted on the same graph in Figure 13. The data of Van Atta and Chen are from 
an atmospheric boundary-layer experiment conducted at sea level. One can see that the 
calculations agree quite well with the experiment suggesting that the shape given by equation 
(1) is a good approximation for the first derivative pdf’s. 

While still on the subject of the first derivatives, we would like to examine the pdf’s 
/«„(«,»;», within the sublayer. The pdf's at the two y + locations within the sublayer 
for both the Reynolds numbers 318 and 395 are plotted in Figure 14. Gaussian pdf is also 
plotted for reference. It can be seen that n, y pdf’s are quite close to Gaussian within the 
sublayer. 

5 Second Derivatives of the Velocities 

We now turn our attention toward the second derivatives of the fluctuating quantities. In 
Figure 15 the pdf’s /«,„(«,**; y, R e — 318) are plotted for the three locations beyond the 
buffer region. It can be seen that all three pdf’s fall on the same curve, indicating that the 
pdf’s outside the viscous sublayer are independent of the Reynolds number. Therefore, a 
single pdf can be calculated by using the data from all of these locations for each of the second 
derivatives. In Figure 16, the pdf’s of the u velocity derivatives u yzx y ^yzy y ^yzz y ^^yyx y ^yyy y ^yxx 
are plotted. On the same plot the pdf’s of the same second derivatives of v and w are also 
plotted. These pdf’s are so close to each other that they are virtually indistinguishable from 
each other. This indicates that similar to the first-derivative pdf’s, the second derivative 
pdf’s are also independent of the y location and all of them fall on the same curve. 

These observations are verified again using the data at R e = 395. The pdf’s, f v , yx {v , yx ; y, /? e )> 
of v, yx plotted in Figure 17 show tlx; independence of the y location. The pdf’s of the same 
second derivative quantities as in Figure 16 at /? f =318 are plotted in Figure 18 at 7? e =395. 
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Ibis again shows the independence of the y location and that the pdf’s fall on the same 
curve irrespective of the quantity. Least square fits to the pdf’s of the second derivatives at 
the two Reynolds numbers are plotted in Figure 19 for comparison. 

It should be mentioned that in general the pdf’s calculated from the data at R c = 395 
exhibit a larger spread along the tails. Most likely this is so because of the relatively smaller 
number (about 100K) of points used in calculating the pdf’s from the data at this Reynolds 
number. 


6 Dissipation Rate 


In this section, we examine the log normality of the kinetic-energy dissipation rate. The 
dissipation rale is defined as 


c 


(ij - 1,2,3), 


( 2 ) 


where repeated indices imply summation and ll e is the Reynolds number of the flow. In 
Figures 20 and 21, the pdf’s of log(e) are plotted for y > y a . On the same plot the Gaussian 
(normal) pdf is plotted for comparison. This clearly demonstrates that at both the Reynolds 
numbers, the dissipation rate is log normal outside the sublayer and the buffer regions. 
However, in the sublayer and the bufTer region the pdf (not shown) of the logarith m of e is 
skewed to the left. 


T Concluding Remarks 

By analyzing the data at the two Reynolds numbers (318 and 395), it is concluded that 
the pdf’s of the fluctuating quantities are independent of y location when outside the bufTer 
layer. The pdf s of u are skewed slightly to the left. The pdf’s of v and w are symmetric 
and appear to be close to Gaussian. 

The pdf’s of the first derivatives are independent of y for locations outside the buffer 
layer. Not only arc the derivative pdf’s independent of y, but also they appear to have the 


same universal shape irrespective of the particular quantity: the shape of a Gaussian core 
and exponential tails. The fact that the calculated pdf’s from a channel flow agree with 
those from a boundary layer experiment suggests a universality of even more general nature 
when sufficiently far away from the boundaries. Further, the pdf’s of u,„ within the viscous 
sublayer seem to be independent of y location and appear close to Gaussian. 

The pdf’s of the second derivatives are also independent of y location and, irrespective 
of the quantity, they all fall on a universal curve. This curve is also made up of a Gaussian 
core and exponential tails, but it is different from the one for the first derivatives. 

The log normality of the dissipation rate for y location outside the viscous sublayer is 
established at both the Reynolds numbers. 

Since all the above conclusions arc verified at two different Reynolds numbers, the next 
question is the Reynolds number independence. This will be examined in a future study. 
The universality of pdf’s in compressible flow will also be investigated. 
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Figure 10. PDFs of the first derivatives for y beyond the viscous sublayer 

for Re=31 8. Symbols are from DNS and the solid line is Gaussian. 










Figure 13. Comparison of DNS data with the experiment of VanAtta and Chen. 

(--) and (...) are the least squares fits to the pdfs of the first derivatives 

from the DNS data at Re=318 and 395, respectively. Symbols indicate 
experiment and the solid line is Gaussian. 
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Figure 14. PDFs of u, y for y within the viscous sublayer. Symbols are from 
the DNS calculations and the solid line is Gaussian. 
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